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Critical exponents in U(l) lattice gauge theory with a monopole term * 
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We investigate critical properties of the phase transition in the four-dimensional compact U(l) lattice gauge 
theory supplemented by a monopole term for values of the monopole coupling A such that the transition is of 
second order. It has been previously shown that at A = 0.9 the critical exponent is already characteristic of a 
second-order transition and that it is different from the one of the Gaussian case. In the present study we perform 
a finite size analysis at A = 1.1 to get information wether the value of this exponent is universal. 
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1. INTRODUCTION 

Our investigations are based on the Wilson ac- 
tion supplemented by a monopole term 

s = p (i-co S e^ x ) + \J2\M P M (i.i) 



p.x 



where M p . x = e p(Tfil ,(Q fl , VjX+a - 



-0|u/,x)/47r and the 
physical flux Q^.x E [— n 7 ir) is related to the 
plaquette angle @u V , x e (-47r,47r) by <d^ x = 



Qfj,v,x + 2im pv . x |2j]. We use periodic boundary 
conditions. 

Studies of energy distributions with a newly de- 
veloped dynamical algorithm have indicated that 
with increasing A the strength of the first or- 
der transition decreases such that the transition 
ultimately gets of second order [HQ]. Recently 
the behavior at larger A has been investigated in 
more detail [|5[ . The unexpected phenomenon has 
been revealed that the phase transition persists 
up to very large values of A, where the transition 
moves to large negative 0. Further it has been 
observed that (within errors) the monopole den- 
sity becomes constant in the second-order region. 
By a finite size analysis it has been shown that at 
A = 0.9 the critical exponent is already character- 
istic of a second-order transition. Moreover, with 
the result v — 0.446(5) obtained, it turned out to 
be different from the exponent h of the Gaussian 
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Similar results have been recently reported || 
for the Wilson action extended by a double charge 
term with coupling 7. There the first order transi- 
tion, which occurs at 7 = 0, weakens with increas- 
ing 7 until it becomes of second order at a tricrit- 
ical point. While for the usual periodic bound- 
ary conditions the second order region starts at 
negative 7, for a spherelike lattice it occurs al- 
ready at 7 = ||. This could be due to inhomo- 
geneities which tend to weaken the transition || . 
The result obtained in Ref. § for 7 = -0.2, -0.5 
is v = 0.365(8) which is smaller than the value 
mentioned above. Thus the two types of actions 
obviously lead to different results. 

In the present study we continue the investiga- 
tions with the action (1.1) because the addition 
of a monopole term is attractive in view of the 
close relation of monopoles to the phase struc- 
ture. The increase of the finite-size effects with A 
observed in || has hinted at a possible increase 
of v. This suggests to check if for larger A the 
same exponent is obtained, which is done here by 
performing a finite size analysis at A = 1.1. 

2. DETERMINATION OF THE TRAN- 
SITION POINT 

In order to facilitate the determination of the 
transition point we have used the topological 
characterization of the phases HQ. It is based 
on the fact that there is an infinite network of 
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Figure 1. Specific heat C as function of (5 for 
A = 0.9 and 1.1 on the 8 4 lattice. 



monopole current lines in the confining phase and 
no such network in the Coulomb phase. On finite 
lattices "infinite" is to be defined in accordance 
with the boundary conditions For the peri- 
odic boundary conditions we are considering here, 
"infinite" is equivalent to "topologically nontriv- 
ial in all directions" . While for loops the topo- 
logical characterization is straightforward, to de- 
termine the topological properties of networks it 
is necessary to perform a more elaborate analysis 
based on homotopy preserving mappings Hfjj. 

Because the probability to find a network which 
is nontrivial in all directions takes values exactly 
1 and in the confining phase and in the Coulomb 
phase, respectively, it is a very advantageous or- 
der parameter. Indeed, with an infinite system 
a single configuration would be enough and on a 
finite lattice still few configurations are sufficient 
to discriminate between the phases. 

Since on finite lattices different order parame- 
ters lead to slightly different critical (3, to deter- 



mine the location of the maximum of the specific 
heat for larger A in an efficient way we first de- 
termine the critical (3 from the topological order 
parameter and then find the maximum of the spe- 
cific heat in an easy second step. 

That finite size effects increase with A is in- 
dicated by the fact that the transition region be- 
comes broader. From Figure 1 it can be seen that 
the width of the peak of the specific heat increases 
and that its height decreases with increasing A. 

3. CRITICAL BEHAVIOR 

In order to obtain information on the nature 
of the phase transition we have investigated the 
finite-size scaling behavior of the maximum of the 
specific heat C max - It is expected to be 

C m ax ~ L (3-2) 

if the phase transition is of first order and 

Cmax ~ L " (3-3) 

if it is of second order, where a is the critical 
exponent of the specific heat and v the critical 
exponent of the correlation length. 

In Figure 2 we present results of simulations for 
C max on lattices with L = 6, 8, 10, 12, for A = 

0. 9 obtained in Ref. || and for A = 1.1 obtained 
here, at the respective values of /3 cr . The fit to 
these data gives the values for — shown in Table 

1 . They are clearly quite far from 4 and thus the 
transition is not of first order. 

Table 1 



A 


0.9 


1.1 


a 
v 


0.485(35) 


0.289(68) 


V 


0.446(5) 


0.466(8) 


&(oo) 


0.4059(5) 


0.1882 (34) 


a 


-1.99(6) 


-2.29 (10) 



The results for v listed in Table 1 are obtained 
using the hyperscaling relation a = 2 — d v. They 
are different from the value \ of the Gaussian 
case. There is a slight increase with A which is 
beyond statistical errors. 
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Figure 2. Maximum of specific heat C max as func- 
tion of lattice size L for A = 0.9 and 1.1 at /3 cr . 



The critical /3 is expected to behave as 

p a {L) = A*(oo) + aL-i (3.4) 

From this relation, using the values of v in Table 
1 and our data for /3 cr (£), the numbers given for 
/3 cr (oo) and a in Table 1 are obtained. 

4. DISCUSSION 

The slight increase of v with A observed could 
indicate a nonuniversal behavior or even hint at 
the possibilty that v ultimately reaches the Gaus- 
sian value. However, finite size analyses on larger 
lattices and for still larger A appear necessary to 
decide these questions. 

Different values for v here and in Ref. || could 
signal different universality classes for the re- 
spective actions. One should, however, also be 
aware of the possibility of obtaining different lim- 
its when approaching from different regions of 
coupling space M. Further the use of different 



operators can lead to different values. Indeed, in 
an investigation of gauge balls with the double 
charge action for 7 = —0.2 the values v ~ 0.35 



and 0.5 have been reported [10 1. 

It is interesting that recently theoretical argu- 
ments have been presented [jllj predicting the set 
of exponents 1/3, 5/11 and 1/2, which within er- 
rors are in agreement with the values observed in 
the simulations in Rcfs. j|||[l(| and here. 

With the non-Gaussian exponents continuum 
limits of four-dimensional pure U(l) gauge the- 
ory, nontrivial and not asymptotically free, ap- 
pears possible. They would be different from the 
usual expectations for QED and would not be 
seen in perturbative approaches. Possibly then 
monopoles survive in the limit. 
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